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We have calculated the effects of structural distortions 
of armchair carbon nanotubes on their electrical transport 
properties. We found that the bending of the nanotubes de- 
creases their transmission function in certain energy ranges 
and leads to an increased electrical resistance. Electronic 
structure calculations show that these energy ranges contain 
localized states with significant a-n hybridization resulting 
from the increased curvature produced by bending. Our cal- 
culations of the contact resistance show that the large contact 
resistances observed for SWNTs are likely due to the weak 
coupling of the NT to the metal in side bonded NT-metal 
configurations. 

Carbon nanotubes (NTs) can be metallic or semicon- 
ducting. They have high mechanical strength and good 
thermal conductivity pL properties that make them po- 
tential building blocks of a new, carbon-based, nanoelec- 
tronic technology fl-|f|. Conduction in defect-free NTs, 
especially at low temperatures, can be ballistic, thus in- 
volving little energy dissipation within the NT Fur- 
thermore, NTs are expected to behave like quasi one- 
dimensional systems (Q1D) with quantized electrical re- 
sistance, which, for metallic armchair nanotubes at low 
bias should be about 6 kf2 (h/Ae 2 ). The experimentally 
observed behavior is, however, quite different. The con- 
tact resistance of single-wall nanotubes (SWNTs) with 
metal electrodes is generally quite high. Furthermore, 
at low temperatures a localization of the wavefunction 
in the nanotube segment contained between the metal 
electrodes is observed that leads to Coulomb blockade 
phenomena Q]. The latter observation suggests that a 
barrier or bad-gap develops along the NT near its contact 
with the metal. In an effort to understand the origin of 
these discrepancies we have used Green's function tech- 
niques to calculate the effect of the modification of the 
NTs by bending on their electronic structure and electric 
transport properties. We also investigated the effects of 
the strength of the NT-metal pad interaction on the value 
of the contact resistance. 

Most discussions on the electronic structure of NTs 
assume perfect cylindrical symmetry. The introduction 
of point defects such as vacancies [|| or disorder ||||] 
has been shown to lead to significant modification of 
their electrical properties. Here we focus on the ef- 
fects of structural (axial) distortions on the transport 



properties of armchair NTs. AFM experiments [|T0[ and 
molecular mechanics simulations jflj have shown that 
the van der Waals forces between NTs and the substrate 
on which they are placed can lead to a significant de- 
formation of their structure. To maximize their adhe- 
sion energy the NTs tend to follow the topography of 
the substrate pO| , pT[ . Thus, for example, NTs bend to 
follow the curvature of the metal electrodes on which 
they are deposited. When the strain due to bending ex- 
ceeds a certain limit, kinks develop in the nanotube struc- 
ture pO| , p^ , ^3| . It is important to understand how these 
NT deformations affect the electrical transport proper- 
ties of the NTs. Could they be responsible for the low 
temperature localization observed ? |jj Early theoretical 
work on this issue was based on tight-binding model in- 
volving only the 7r-electrons of the NTs and accounted 
for the electronic structure changes induced by bending 
through the changes in 7r-orbital overlap at neighbor- 
ing sites. This study concluded that bending distortions 
would have a negligible effect on the electrical proper- 
ties of the NTs [^i). The applicability of this approach 
is limited to weak distortions. Experiments, however, 
show that strong deformations and kink formation are 
common. Under such conditions, bending-induccd o-~k 
mixing, which was not considered before, becomes very 
important in strongly bent NTs In this work, the 

NT electronic structure is computed using the extended 
Hiickel method (EHM) jljj that includes both s and p va- 
lence electrons. We have previously jl6| shown that EHM 
calculations on an armchair (6, 6) NT model (96 A long) 
reproduce the electronic properties obtained with more 
sophisticated ab-initio and band structure computations 
on NTs. The approach we used in the computation of 
the electrical properties is similar to that of Datta et al. 

The conductance through a molecule or an NT cannot 
be easily computed; even if the electronic structure of the 
free molecule or NT is known, the effect of the contacts 
on it can be substantial |^0| and needs to be taken into 
account. Typically, there will be two (or more) leads con- 
nected to the NT. We model the measurement system as 
shown in Figure la. The leads are macroscopic gold pads 
that are coupled to the ends of the NT through matrix 
elements between the Au surface atoms and the end car- 
bon atoms of the NT. In most experiments to date the 
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NTs are laid on top of the metal pads. As we discussed 
above, the NTs then tend to bend around the pads. Such 
bending deformations are modelled in our calculations by 
introducing a single bend placed at the center of the tube. 

The electrical transport properties of a system can 
be described in terms of the retarded Green's function 
p8| , pH . To evaluate the conductance of the NT we need 
to compute the transmission function, T(E), from one 
contact to the other. This can be done following the 
Landauer-Biittiker formalism as described in |18|. The 



key element of this approach lies in the treatment of the 
infinite leads which are here described by self-energies. 
We can write the Green's function in the form of block 
matrices separating explicitly the molecular Hamilto- 
nian. After some simplification we obtain: 



Gnt = [ESnt ~ Hnt — Si — £2] 



(1) 



where Snt and Hnt are the overlap and the Hamilto- 
nian matrices, respectively, and £1,2 are self-energy terms 
that describe the effect of the leads. They have the form 
t} gjTj with gi the Green's function of the individual leads 
fl9|,P2l and Tj is a matrix describing the interaction be- 
tween the NT and the leads. The Hamiltonian and over- 
lap matrices are determined using EHM for the system 
Gold-NT-Gold. The transmission function, T(E), that is 
obtained from this Green's function is given by pSfl : 

T(E) = T 21 = Tt^GntTxG^]. (2) 

In this formula, the matrices have the form: 

r 1)2 = t(Ei, 2 - sj j2 ). (3) 

The summation over all conduction channels in the 
molecule allows the evaluation of the resistance at the 
Fermi energy, R = h/ (2e 2 T(Ep)). Transport in the pres- 
ence of an applied potential is also computed. The dif- 
ferential conductance is computed in this case using the 
approximation jlq]: 



dl 2e 2 



(1-1)%)] 



(4) 



with i] G [0, 1] describing the equilibration of the nan- 
otube energy levels with respect to the reservoirs [fL9| . 
As a reference, we use the Ef obtained from EHM for 
individual nanotubes as the zero of energy. The NT 
model used in our calculations is a (6,6) carbon nan- 
otube segment containing 948 carbon atoms. The bond 
distance between carbon atoms in non-deformed regions 
of the NT is fixed to that in graphite 1.42 A, leading to a 
tube length of 96 A. The building of deformed NTs using 
molecular mechanics minimization schemes |23|p4[ has 
been described in detail elsewhere . The structures of 
the bent NTs are shown in Figure lb. The metallic con- 
tacts consist each of 22 gold atoms in a (111) crystalline 



arrangement. The height of the NT over the gold layer 
is 1.0 A, where the Au-C bond distances vary from 1.1 
to 1.6 A. 
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FIG. 1. (a) Schematic of the model used in computation 
involving a finite length carbon nanotube contained between 
infinite gold pads, (b) Structures of the isolated bent nan- 
otubes used in the computation. The indices refer to the 
bending angles (in degrees). 

In Figure 2 we present the computed transmission 
function T(E) for the bent tubes (note that T(E) rep- 
resents the sum of the transmission probabilities over all 
contributing NT conduction channels). The upper- right 
of Figure 2 shows the raw transmission results obtained 
for the straight NT. The fast oscillations of T(E) are due 
to the discrete energy levels of the finite segment of the 
carbon nanotubes used. For clarity, we will use smoothed 
curves in the description of the results. At Ep, T(E) is 
about 1.2, leading to a resistance (ss 11 kil ) higher than 
expected for ballistic transport (« 6 kil for T(E) — 2.0). 
This reduction in transmission is due to the contribu- 
tion from the contact resistance. The increasing T(E) 
at higher binding energies is due to the opening of new 
conduction channels. The asymmmetry in the transmis- 
sion T(E) is a function of the NT-pad coupling (C-Au 
distance). A longer NT-Au distance increases the T(E) 
above Ep, while it decreases it below Ep, and vise versa. 
Since the NT-pad geometry is kept fixed in all computa- 
tions, this behavior does not influence the effects induced 
by NT bending. 

According to our calculation the contact resistance at 
Ep is only about 5 kf2, much smaller than the rj 1 Mf2 
resistance typically observed in experiments on single- 
wall NTs [g||. The dependence of T(E) and contact 
resistance at Ep on the Au-NT distance is shown in 
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the upper-left of Figure 2. We see that T(Ep) remains 
nearly constant between 1-2 A, then decreases exponen- 
tially. For distances appropriate for van der Waals bond- 
ing (> 3 A) the contact resistance is already in the Mf2 
range. The above findings suggest that the high NT con- 
tact resistance observed experimentally may, in addition 
to experimental factors such as the presence of asperities 
at the metal-NT interface, be due to the topology of the 
contact. In most experiments, the NT is laid on top of 
the metal pad. The NT is at nearly the van der Waals 
distance away from the metal surface, and given that 
transport in the NT involves high fc-states which decay 
rapidly perpendicular to the tube axis, the coupling be- 
tween NT and metal is expected to be weak p5fl . Direct 
chemical bonding between metal and the NT, or interac- 
tion of the metal with the NT cap states 1 26 1 should lead 
to stronger coupling. In this respect, it has been found 
pjj that high energy electron irradiation of the contacts 
leads to a drastic reduction of the resistance. Since the 
irradiation is capable of breaking NT C-C bonds it may 
be that the resulting dangling bonds lead to a stronger 
metal-NT coupling. 
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FIG. 2. Transmission function T(E) for the different bent 
nanotubes. The upper-right inset gives the raw data obtained 
from the computation for the straight NT along with the 
smoothed curve. The upper-left shows the variation of the 
transmission function T(Ef) with the nanotube-pad distance 
(R(NT-pad)). In this last figure, one gold-NT distance is fixed 
at 1.0 A while the other is varied. 

The strongest modification of T(E) as a result of bend- 
ing is observed at around E=-0.6 eV where a transmis- 
sion dip appears. This dip is strongest in the 60° bent 
NT. Furthermore, its transmission function at higher 



binding energies (BE) is lower than those of the 0°-45° 
bent NTs, indicating that the transmission of higher con- 
duction channels is also decreased. The nature of the dip 
at about -0.6 eV can be understood by examining the 
local density-of-states (LDOS) of bent tubes shown in 
Figure 3 |15|. A change (increase) in the LDOS is seen 
in the same energy region (0.5-0.8 eV below Ep) as the 
transmission dip. This change is essentially localized in 
the vicinity of the deformed region. The new states result 
from the mixing of a and 7r levels, have a more localized 
character than pure n states leading to a reduction of 
T(E). As Figure 3 shows, the change in transmission 
with bending angle is not gradual; the transmission of 
the 30° and 45° models is only slightly different from 
that of the straight tube. Apparently, large changes in 
DOS and T(E) require the formation of a kink in the NT 
structure, as is the case in the 60° and 90° bent NTs. 
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FIG. 3. Variation of the LDOS near the Fermi level for 
several bent (6, 6) armchair nanotubes. (resolution = 0.05 
eV). Indices give the relative position of the section in the 
nanotube structure (1: boundary, 40: middle of the nanotube) 

Once the transmission function is computed, the deter- 
mination of the differential conductance and resistance is 
straigtforward. Figure 4 shows the results for two ex- 
treme cases of equilibration of the Fermi levels. The 
first is when r\ = (Figure 4a), and the symmetric case 
i] = 0.5 (Figure 4b). When ?y=0.0, the Fermi level of 
the NT follows exactly the applied voltage on one gold 
pad and the conductance spectrum is directly propor- 
tional to T(E). As expected, there is no large difference 
between the 0°, 30° and 45° models, while the 60 and 
90° models show the dip structure at around 0.6 V. The 
non-linear resistance (NLR) spectra show clearly a sharp 
increase by almost an order of magnitude at 0.6 V. These 
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features are also observed when 77=0.5, where the Fermi 
level of the NT is floating at half the voltage applied be- 
tween the two gold pads. The dip at around 1.2 V in 
the conductance spectra is now broader, and the NLR of 
the 60° bent tube increases by about a factor of 4 from 
the computed resistance of the straight tube. These re- 
sults suggest that there exists a critical bending angle 
(between 45° and 60°) above which the conduction in 
armchair carbon nanotubes is drastically altered |2£| . 



that NTs end-bonded to metal pads will have contact re- 
sistances of only a few kf2. Such low contact resistances 
will greatly improve the performance of NT-based devices 
and unmask the Q1D transport properties of NTs. 
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FIG. 4. Differential conductance (top) and resistance (bot- 
tom) of bent tubes for two extreme cases where (a) 77=0.0 and 
(b) 77=0.5 

In conclusion, we have calculated the effects of struc- 
tural distortions of armchair carbon nanotubes on their 
electrical transport properties. We found that bending of 
the nanotubes decreases their transmission function and 
leads to an increased electrical resistance. The effect is 
particularly strong at bending angles higher than 45° de- 
grees when the strain is strong enough to lead to kinks 
in the nanotube structure. The electronic structure cal- 
culations show that the reduction in T(E) is correlated 
with the presence at the same energy of localized states 
with significant a-n hybridization due to the increased 
curvature produced by bending. Resistance peaks near 
Ep are the likely cause for the experimentally observed 
low temperature localization in carbon NTs bent over 
metal electrodes (7). Our calculations of the resistance 
(including the contact resistance) of a perfect NT give a 
value close to ft./2.4e 2 instead of h/Ae 2 . This increase in 
resistance is solely due to the finite transmission of the 
contacts. The much larger contact resistances observed 
in many experiments on SWNTs are likely due to the 
weaker coupling of the NT to the metal when the NT is 
simply placed on top of the metal electrodes. We predict 
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